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Is soft physics entropy driven? 
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Abstract. The soft physics, pt < 2 GeV/c, observables at both RHIC and the SPS have now been mapped 
out in quite specific detail. From these results there is mounting evidence that this regime is primarily 
driven by the multiplicity per unit rapidity, dNc^/dr?. This suggests that the entropy of the system alone is 
the underlying driving force for many of the global observables measured in heavy-ion collisions. That this 
is the case and there is an apparent independence on collision energy is surprising. I present the evidence 
for this multiplicity scaling and use it to make some extremely naive predictions for the soft sector results 
at the LIfC. 
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1 Introduction 

Back in 1965, Hagedorn postulated that the hadronic mass 
spectrum grows exponentially with mass pp. Since 1965, 
more than 3200 additional resonance states have been 
identified and these results seem to confirm Hagedorn's 
hypothesis. The exponential growth means that adding 
more energy to a system merely produces more particles, 
resulting in a limiting temperature, Tq ~ 160 MeV for 
a hadron gas. Interestingly this value of Tq is very close 
to the limiting chemical freeze-out temperature, Tch, cal- 
culated via statistical production models, from data as a 
function of ^/s^, (see Fig.^ 2^). We know that entropy 
is related to temperature and energy via TAS=AE. If the 
calculations are correct and chemical freeze-out tempera- 
tures reach a maximal value then, as the collision energy 
increases, the entropy of the system must also increase. 

Even before Hagedorn's work Landau and Fermi re- 
lated pion production in high energy collisions to the en- 
tropy produced TT. I will try and reproduce the main 
points of their arguments here. The energy density, e, 
available for particle creation in a collision is 



V ~ 2mArVb 



(1) 



where: Vq is the volume and m^r is the mass of the nu- 
cleon. If you assume that the entropy is produced early 
in the collision, that the source is thermalized and that 
it expands adiabatically, the entropy, S, can be related to 
e via the equation of state (EoS). Thus, entropy and the 
produced particle multiplicity are correlated. 

A simple example of such a correlation can be derived 
using the EoS of a system of massless pions which emits 
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Fig. 1. The calculated Tch as a function of collision energy, 
for e"'"-e~ (squares), p-p (open triangles) p-p (circles) and A- A 
(closed triangles). Data from ^3- 



as a blackbody. In this case, p = e/3 and e = T"*. Thus 

s = S/V 



Ts. 



+ P 
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(2) 
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S ~ V{{y/SNN - 2mN)y/SNN)^^^ 
■1/4 



\/snn 



(5) 
(6) 



The ^s„„ dependence of the entropy in Eqn. is the 
same as that derived by Fermi ^ and S/V is called the 
Fermi energy, F. If these arguments hold, a plot of the 
measured entropy density, or mean multiplicity per par- 
ticipant nucleon, versus F, calculated from the right-hand 
side of Eqn. |^ should result in a linear correlation. Such 
a correlation can be seen in Fig. [3 taken from j^. In 
this plot, although only the mean pion multiplicities are 
used, this is a good approximation to the total multiplic- 
ity. It can be seen that the multiplicity per participant of 
the p-p collision data, open circles, is indeed proportional 
to F. The A- A data seem to depart from this curve at 
the top AGS energy. This apparent increase in entropy 
in higher energy collisions was first observed by 1995 and 
interpreted as evidence of a transition to quark degrees of 
freedom resulting in a different EoS ^ . 
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Fig. 2. Entropy per participant versus the Fermi energy vari- 
able (only pion multiplicities were used) for various collision 
energies 0. 



Since there appears to be evidence that particle multi- 
plicity and entropy are correlated over a range of collision 
energies, it is interesting to see if entropy, e.g. multiplicity, 
can be used to determine unique scalings for other bulk 
physics measurements. The rest of this paper discusses 
femtoscopic measurements, elliptic flow and strangeness 
yields and uses entropy scalings to make naive predictions 
for these observables at the LHC. 



2 Femtoscopy 



Femtoscopic measurements examine space-momentum cor- 
relations and allow, in a model dependent fashion, for 
the extraction of source sizes. For more details on these 
techniques see elsewhere 8j. Detailed measurements have 
been made, especially looking into identical pion correla- 
tions. While the radii hover around 5 fm, when the re- 
sults are plotted as a function of ^/s^, see Fig. left 
panel no clear trends appear. However, if, as in Fig. |31 
right panel, the same measurements are now plotted as 
a function of (dNch/dr/)^/'^ a clear scaling is observed for 
Hiong and Kside- The applicability to Rom* is less clear. 
This lack of scaling for Rout might have an explanation 
via the fact that this HBT variable consists of a convolu- 
tion of terms containing length and time scales. R;ong and 
Rside are, on the other hand, pure length measurements. 
The deviations in Kout might therefore be signatures of 
differing timescales for the collisions. It should be noted 
that while the right panel of Fig. is for < kx > — 400 
MeV and 390 MeV for RHIC and SPS respectively. As 
good a scaling is observed if other < kx > ranges are se- 
lected. The non-monotonic behavior seen in the left panel 
of Fig. 121 for y/Sj^„ < 5 GeV makes it obvious that the 



dNc;i/d?7 scaling breaks down at very low collision ener- 
gies. It has been pointed out, by the CERES collaboration, 
that this is probably due to the dominance of baryons at 
lower 
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Fig. 3. Left panel) World dataset of published HBT radii 
from central Au-|-Au (Pb-|-Pb) collisions versus collision en- 
ergy. Compilation from [Jj. Right panel) The same measure- 
ments as in the left panel but as a function of [dNch/dr])^^^ 
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3 Elliptic Flow 

Peripheral and low energy collisions are likely to produce 
systems with incomplete thermalization. Since re-scattering 
of the particles is rare, in this low density limit regime, 
little change occurs, on average, to the initial momentum 
distributions. The measured elliptic flow, V2, is therefore 
proportional to the initial state eccentricity, e, as deter- 
mined via Glauber calculations, and the space density of 
the initial particles. The latter determining the scattering 
probability which must be finite so that the initial spatial 
anisotropy can be converted into a measurable momentum 
anisotropy. Thus, 



V2 (X 



1 dN 
Area dy 



(7) 
(8) 



where Ra, and Rj, are the major and minor axes of 
the overlap ellipse of the colliding nuclei in the transverse 
plane. From Eqn.|Slit can be seen that in the most central 
A-A collisions e, and hence V2, tends towards zero. In the 
Cu-Cu data however there remains a significant measured 
V2 component in the most central events. Fig. 01 
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Fig. 4. Elliptic flow, V2, measured near mid-rapidity as a func- 
tion of Npart for Cu-Cu (open symbols) and Au-Au (closed 



symbols) collisions at 



= 200 GeV. Only statistical er- 



rors are shown, taken from 



Also when, for a given number of participants, Npa^t, 
V2/e measurements in Cu-Cu are contrasted to similarly 
scaled Au-Au results, they are not compatible. This sug- 
gests that the low density limit concept is not applicable 
to the data Since the Cu-Cu results are above those 
of the Au-Au at a given Npa^t , other conclusions could be 
drawn from these data. These are that either the Cu-Cu 
system is more efficient at translating spatial anisotropics 
into V2 signals or there are significant non-flow effects in 



the Cu-Cu data that create these large momentum asym- 
metries. An alternative explanation can however be pos- 
tulated by examining the initial state eccentricity calcula- 
tion. As stated above this is traditionally defined as an av- 
erage eccentricity calculated from the participant nucleon 
distributions relative to the reaction plane for each cen- 
trality class via Glauber calculations. This method leads 
to two possible errors: firstly event-by-event fiuctuations 
are averaged to zero and secondly, while the minor axis 
of the participant distribution is along the impact param- 
eter vector on average for any given event this will not 
generally be the case. In events with a small number of 
participants these two effects can cause a significant error 
in the eccentricity calculations The PHOBOS col- 
laboration have therefore proposed a new technique for 
calculating the eccentricity, called tpart , which takes these 
effects into account ^^I- When the Y2/^part ratio is then 
plotted versus the mid-rapidity area density of produced 
particles, Fig.jS] good agreement is now observed between 
the two collision systems. Also shown in this plot are re- 
sults from yfs^^ = 17 and 4 GeV and a common scaling 
seems to emerge. This suggests that V2 depends solely on 
the initial area density achieved in the collision over two 
orders of magnitude in collision energy. Note that a linear 
dependence on dN/dy or entropy has emerged. 



03 - 



0.1 - 



1 

PHOBOS prelrminary 

• 200 GeV, Au, trks 

200 GeV, Au, hits 

130 GeV, Au, hits 

200 GeV, Cu, t*s 

200 GeV, Cu, hits 

62,4 GeV, Cu, hits 



130 GeV, Star 
17 GeV, Na49 
4 GeV, £877 




10 



20 



30 



1/<S)<dN^^/dy> [fm^] 



Fig. 5. V2 per average participant eccentricity as a function 
of mid-rapidity area density. Data are from RHIC Au-Au col- 



lisions at 



=130 and 200 GeV and Cu-Cu coUisions at 



./s-^ = 62.4 and 200 GeV as weU as SPS 
and AGS ,/ 



= 17 GeV 



= 4 GeV collisions [TT]- 



4 Strange particle production 

A lack of phase space in p-p collisions is predicted to result 
in strangeness production being suppressed . The mag- 
nitude of this suppression decreases with increasing ^s„„ 
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and with the volume of the medium produced. Since the 
colhsion volume is linearly proportional to Npart, it is also 
predicted that the strangeness production volume follows 
such a linear scaling. The calculations also show that the 
p-p suppression increases with the particle's strangeness 
content. When the system becomes sufficiently large, the 
phase space restrictions are removed and strangeness can 
be produced freely. For these types of collisions, the ratio 
of strangeness yield to Npart is therefore a constant as a 
function of centrality. 
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Fig. 6. Yield per participant relative to the yield per partici- 
pant in p-p collisions as a function of Npart- Solid symbols are 
for Au-Au collisions at ^/s^ — 200 GeV, and the open sym- 
bols are measured by NA57 from Pb-Pb at ^/s'^ =17.3 GeV. 
The boxes represent the combined statistical and systematical 
uncertainties in the p-p and p-Be data. The error bars on the 
data points represent those from the heavy-ion measurement. 
The boxes on the right axes mark the predictions from a model 
using a Grand Canonical formalism described in |13|. 



It can been seen in Fig. that in Au-Au collisions at 
~ = 200 GeV, although the predicted strangeness or- 



dering of the p-p suppression is observed the yields per 
participant are not constant. Further, the magnitude of 
the suppression for a given strange baryon is the same in 
= 200 GeV collisions as at the SPS in ^/sTT = 17.3 



= ^UU UeV collisions as at tne Sl^S m ^Sj^ 
GeV collisions [TmT3] . Both of these results are counter to 
the predictions. Calculations using statistical models ap- 
plied to the STAR data indicate that 7^, the strangeness 
saturation factor, is unity in central collisions, e.g. that 
there is no suppression in the most central data IQ^ . This 
suggests that the relevant volume for strangeness pro- 
duction is not linearly proportional to Npa^t and hence 
not purely related to the initial collision overlap volume. 
Therefore some other scaling must be sought. 
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Fig. 7. Yield per dNch /d?) as a function of dN^h /dr] for strange 
(anti)baryons. Solid symbols are for Au-Au collisions at y^s^^ 
— 200 GeV, and the open symbols are measured by NA57 from 
Pb-Pb at ^s;;^ =17.3 GeV [TK]. 



Figure [3 shows the yield per dNch/drj for A, A, ^ , 
E:+, and f2-+f2+ for RHIC and the SPS Pb-Pb collisions 
using data taken from |14II15| . When plotted in this way 
the yield/dNch/d?/ appears constant for the more central 
data suggesting that an entropy measure is more closely 
correlated to the strangeness production volume than the 
number of participants. 



5 LHC predictions 

As shown in Fig. [3 it appears that the EoS for heavy-ion 
collisions is different to that in p-p. We also expect it to 
be much more complicated than that derived earlier for 
a gas of massless pions. There can still be a functional 
relationship between entropy and ^s„„. In Fig. |S1 taken 
from |17j , measurements of the mid-rapidity dNch /dt] per 
participant pair are shown as a function of collision en- 
ergy, note the log scale on the abscissa. This plot clearly 
shows a linear scaling of Nch/drj per participant pair with 
log(^s„„) again over several orders of magnitude. This 
dramatic scaling allows an estimation of the change in en- 
tropy as a function of beam energy to be made. 

Making a leap of faith, one can naively extend this scal- 
ing to the LHC energy regime of 5.5 TeV, as performed in 

and shown by the sohd line in Fig.|Sl Such an extrap- 
olation leads to the expectation that mid-rapidity multi- 
plicities will be as little as 60% higher than they are at 
RHIC, e.g. the mid-rapidity yield will be of order 1200 
charged particles. Armed with this "knowledge" of the 
expected LHC multiplicity/entropy one can continue to 
make simple extrapolations. I use the scalings with en- 
tropy for the other bulk physics measurements discussed 
above. Thus predictions for the HBT radii, elliptic flow 
and strangeness yields are obtained. 
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Fig. 8. The multiplicity density per participant pair measured 
in Au-Au and Pb-Pb coUisions at AGS, SPS, and RHIC. Taken 
from |17|. and naively extrapolated to LHC energies in 



The femtoscopic radii were shown to depend primarily 
on the cube root of the mid-rapidity event multipHcity. 
Thus the combination of Figures 01 and |H1 suggests that 
HBT measurements in central collisions at LHC will be 
~ 17% higher than at RHIC ((1.6)^/3 ^ ;l.17). Precise ex- 
pectations for Rout at the LHC are, admittedly, less cer- 
tain, as the former does not scale exactly with multiplicity 
(Fig.©. 

To make a prediction for the value of elliptic flow in 
the most central collisions one needs a calculation of the 
mid-rapidity transverse area. For Pb-Pb collisions Glauber 
calculations show an initial RMS transverse radius of ~2.8 
fm and hence a mean particle area density of ~ 50 fm~^. 
Using a straight line extrapolation from the data shown 
in Fig. |S1 one would predict V2/epart = 0.35. However, it 
should be remembered that the proportionaUty of ellip- 
tic flow to the mean transverse area density is only pre- 
dicted to hold in the low particle density limit. It is hard 
to believe that central Au-Au or Pb-Pb collisions produce 
low area densities and there is much evidence that the 
re-scattering rate is high at RHIC, for instance the large 
amount of radial flow [181119) and the significant regen- 
eration of resonances after chemical freeze-out Tak- 
ing these, and other, RHIC results into consideration it is 
more likely that central RHIC/LHC collisions are hydro- 
dynamical in behavior. If this is the case, V2/e will be a 
constant, assuming the same EoS holds at RHIC and the 
LHC. If one takes a second look at Fig. 011 is possible that 
the higher centrality Au-Au data points arc pulling away 
from a straight line fit towards a plateau. Therefore in- 
stead of LHC collisions increasing further the elliptic flow 
to eccentricity ratio is likely to stay constant at V2/e = 
0.2. 



Although the strangeness yields cannot be directly ex- 
tracted from Fig. |7| an estimate can still be made for the 
LHC regime. From Fig. it can be seen that whereas 
the A/ A ratio at the SPS is ~0.1, at RHIC the mid- 
rapidity anti-baryon to baryon ratio is very close to unity. 
It can also be seen that the SPS strangeness yields per 
Nc/i/d?7 for the (anti)yl and (anti)H frame those measured 
at RHIC. These results suggest therefore that the LHC 
net-baryon number will be essentially zero and that the 
yield per charged particle at mid-rapidity will be the same 
as that measured at ^/s^ = 200 GeV. Given that the en- 
tropy extrapolations predict an ~60% increase in charged 
particle production one can therefore expect a similar in- 
crease in the strange baryon yields when going from RHIC 
to the LHC. Thus, dN^/dy = dN^/dy = ~ 20 - 30, 
dNH-/dy = dNe+/dy = ~ 4 - 6 and dN,,-+,^+/dy = 
~ 0.5-1 in the most central events. 



6 Summary 

In summary, from simple arguments it can be shown that 
produced particle multiplicities are strongly correlated to 
the entropy of the system and that a limiting temperature 
is reached if the matter being probed is hadronic. Further, 
when heavy-ions are collided at high energies, the trend of 
this entropy measure versus the Fermi energy is different 
to that in p-p and low ^Js^^ heavy-ion collisions. 

Evidence is mounting that many observables from the 
soft-sector are determined primarily by total multiplicity, 
i.e. entropy. The HBT radii, properly-scaled elliptic flow, 
and strangeness yields, all appear to show universal mul- 
tiplicity scalings. This implies that bulk observables are 
dominated by entropy-driven factors. There is also much 
evidence that these high-energy collisions produce a new 
phase of matter that consists of interacting quarks and 
gluons. That this highly non-trivial source apparently pro- 
duces signatures, in the soft momentum sector, with de- 
pendencies only on multiplicity and not reaction energy is 
very intriguing. 

These multiplicity scalings can then be used to make 
naive predictions for some soft sector results at the LHC. 
However, anything can, of course occur when this new 
regime is probed as novel physics could be probed. How- 
ever, any significant deviation from these apparent univer- 
sal multiplicity scalings would be extremely, interesting. 
Whatever the outcome the first results from CERN, now 
only just over a year away, are eagerly awaited. 
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